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Abstract 

It has been suggested that the peak in the specific heat observed numerically for random surface 
actions with extrinsic curvature on dynamical lattices might be the result of a low mass bound state 
in an asymptotically free theory, rather than the signal for a real phase transition. The 0(3) model 
on a fixed lattice displays just such behaviour, but in general transitions appear to be weakened when 
the models concerned are put on dynamical lattices (ie coupled to 2d quantum gravity). We have 
therefore performed simulations of the 0(3) model on dynamical ifi 3 graphs to see if there is still a 
peak in the specific heat and compared the results with those on fixed lattices. 
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There has been a considerable amount of numerical work devoted to exploring the issue of whether a 
non-trivial continuum limit exists for actions of the form 



(1) 



<ij> Ai,Aj 



on dynamically triangulated random surfaces where the ni are the unit normals on neighbouring triangles. 
The theory apparently has a low A crumpled phase and a large A smooth phase similar to those displayed 
by similar actions on fixed triangulated surfaces jl| . The dynamical lattice simulations are of interest for 
constructing well-defined lattice versions of string theory [0, |l 0, ||. The initial simulations in |2j saw 
what appeared to be a second order transition on small lattices, but later work || ||, [| with larger lattices 
and better statistics has shown that the transition, if it exists, is most probably higher order. Indeed, it 
was observed in Q] that the data was not inconsistent with a crossover as the correlation length reached 
the size of the surfaces simulated or with a low mass bound state in an asymptotically free theory. The 
further simulations of Q tend to exclude the first hypothesis, but the second remains a possibility. 

The original one-loop calculation of Polyakov || for the random surface action with extrinsic curvature 
suggested that the extrinsic curvature coupling was asymptotically free and hence there was no phase 
transition. The second (extrinsic curvature) part of the action in equ.(|l|) is, in effect, a sigma model 
living on the dynamically triangulated surface subject to the constraint that the n's be normals to the 
surface. If one forgets these constraints for a surface embedded in three dimensions, which is where most 
of the numerical simulations have been done, it is a 2d 0(3) model 



with \rii\ 2 = 1, which is known to be asymptotically free. The constraints on the normals in Q modified 
the numerical coefficient in the one-loop beta function by comparison with the 0(3) model but did not 
change its sign. A more recent calculation by Polyakov for a surface embedded in four dimensions with 
an action incorporating a topological term has in fact shown that the constraints on the normal vectors 
are softened by renormalization and may thus, in this case, be irrelevant anyway j7|. Intriguingly, a 
calculation of the renormalization flow of the coupling has recently been carried out for an "standard" 
O(N) model coupled to 2d gravit y § where it was found that the result without gravity was renormalized 
by a factor (2/Qa), where Q = y/(26 - N)/3 and a = (V26 - N- V2~N)/Vl2. The rescaling becomes 
complex for the region N > 2 (ie C > 1) where the model without gravity is asymptotically free, which is 
yet another manifestation of the C — 1 barrier in continuum 2d quantum gravity. It is thus unclear what 
to expect from a simulation of the 0(3) model on dynamical lattices from these analytical arguments. 

It was suggested in || that the 2d 0(3) model might serve as a qualitative model for the behaviour of 
discretized random surface actions with extrinsic curvature at the crumpling transition. For a fixed lattice 
2d 0(3) model one observes a peak in the specific heat that increases with system size for small lattices 
before saturating even though the model is known to be asymptotically free. The peak is explained by 
the excitation of a low mass state, the a particle J^, [w| , which contributes when its mass is comparable to 
the inverse correlation length of the lattice 0(3) model. A similar peak is observed in SU(2) lattice gauge 
theory, and it appears in both these cases that there are complex zeroes lurking near the real axis which 
do not converge to it in the infinite volume limit (T^j . In general |l2| transitions appear to be weakened 
when models are transcribed from fixed to dynamical lattices so it is a priori possible that the already 
weak peak in the specific heat of the 2d 0(3) model might get washed out completely when the model is 
put on a dynamical lattice thus invalidating the qualitative parallels that have been drawn between the 
fixed lattice 2d 0(3) model and the dynamically triangulated random surface action of equ.(Q). 

In this paper we conduct a simulation of the 2d 0(3) model on dynamical <fi 3 graphs, which is therefore 
one step closer to the simulations of the surface action in cqu.([l]) than the fixed lattice 0(3) model 
simulations. In our case we are neglecting "only" the gaussian term and the constraints on the vectors 
fii. The simulations can also be regarded as a test of whether anything catastrophic happens to the lattice 
0(N) model for N > 2, as the results of |f| suggest. We choose to employ a simple Metropolis update 
for the spins, which at first sight might appear rather perverse in view of the excellent cluster update 
methods that are now available |l3|. There are two reasons for this choice: Firstly, the simulations of 
equ. (|l|) which we are attempting to approximate employed local algorithms as cluster updates appear to 
be very difficult, if not impossible, to implement for actions containing extrinsic curvature terms especially 
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on dynamically triangulated surfaces; Secondly, a simulation of the 0(3) model on dynamical 3 graphs 
using the Wolff cluster algorithm revealed that there are equilibration problems between the mesh and 
the spins on dynamical triangulations that first appear at quite modest j3 values. This is peculiar to 
dynamical lattices, as a simulation of the 0(3) model on a fixed (f> 3 graph gave identical results for the 
Metropolis and Wolff algorithms. Rather counterintuitively the solution to the problem when using the 
Wolff algorithm appears to be the use of less flip moves per Wolff sweep as (3 increases. 

For the measurements reported here we simulate the 0(3) model on dynamical 3 graphs of spherical 
topology with a fixed number N of points and no tadpole or self-energy insertions which would correspond 
to degenerate triangulations on the direct lattice. We simulate graphs with N from 100 to 10000 points for 
a range of j3 values from 0.05 to 10.0. For each data point we carried out 10,000 metropolis equilibration 
sweeps, followed by 50,000 metropolis measurement sweeps. We used the usual rule of thumb of carrying 
out N flip moves per metropolis sweep, which was sufficient to ensure equilibration between the lattice 
and spin model. We measured the usual spin model properties, the energy 



and susceptibility. The metropolis acceptance was also measured during the course of the simulation and 
the proposed moves adjusted so that it remained in the region of fifty percent as (3 was varied. For the 
lattice itself we measure AL and AF which relate to the acceptance of the flip moves on the graphs. A 
flip can be forbidden either from constraints arising from the graph (ie no tadpoles and no self-energy 
bubbles) or from the energy change in the spin model induced by the reconnection of the vertices. AL 
measures the fraction of randomly selected links which pass the first test and could be flipped according 
to the graph constraints and AF measures the fraction of the links satisfying the graph constraints that 
actually are flipped, ie pass the Metropolis test using the 0(3) model energy change. We also measure 
the fraction of rings of length 3, Pi?3, which serves as an indicator of the local curvature distribution in 
the </> 3 graph. 

Remembering that each point on the lattice is trivalent and taking into account our summation 
conventions we would expect the energy to approach 4/3 for asymptotically large (3. We find that this 
is, indeed, the case for the metropolis algorithm for both fixed and dynamical 4> 3 graphs, although the 
approach is rather slower on dynamical graphs and does not reach E ~ 1.33 until j3 ~ O(10). The Wolff 
algorithm with the number of flips per measurement sweep taken as N fails to reach this asymptotic 
value completely. We are of course principally interested in the behaviour of the specific heat curve on 
dynamical lattices and this is shown in Fig.l for (3 up to 4.0. It is clear from the graph that the specific 
heat looks very similar to that of the fixed lattice 0(3) model. The peak in the specific heat thus still 
mimics rather well the behaviour of the dynamically triangulated random surface action in equ.(l). The 
peak grows for relatively small numbers of points but saturates with larger N - the points for N = 10000 in 
Fig.l are rather poorly equilibrated and are probably best discarded in making comparisons. In addition 
to showing that the 0(3) model on dynamical lattices still produces similar qualitative behaviour to the 
random surface action the simulation provides another example of numerical work venturing into regions 
where continuum analytical calculations suggest trouble (in this case the complex scaling found in Q for 
N > 2) and finding nothing untoward. Simulations and series extrapolation of multiple Ising models on 
dynamical lattices which also have C > 1, at least naively, also fail to show any pathologies [|l4|, whereas 
continuum calculations produce complex critical exponents fl5| . 

The behaviour of the dynamical cf> 3 graphs is shown in Fig. 2. where AF, AL and PR3 are plotted 
for N = 5000. There is a drop in the flip acceptance AL that is associated with the graph constraints at 
the same position as a peak in the probability of rings of length three, PR3, at a j3 value slightly below 
the peak in the specific heat. There is also a dip in the actual flip acceptance AF that is slightly shifted 
from this. These characteristics are again surprisingly similar to those observed in earlier random surface 
simulations, where similar peaks and dips are also manifest. This behaviour is not totally universal for 
all models on dynamical lattices, as AF and AL look rather different for the XY and SOS models |Q, 
though Potts and multiple Potts models on dynamical lattices which have genuine phase transitions give 
results similar to those seen here for a model which presumably does not. 
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In summary, we have added to the repertoire of models that have been simulated on dynamically 
triangulated lattices or their duals by simulating a model that is asymptotically free on a fixed lattice. 
We found that the 0(3) model still possessed a peak in its specific heat that eventually stopped growing for 
large N and displayed very similar lattice characteristics to random surface actions too. These results, 
like the fixed lattice 0(3) results in ||, show that it is not impossible that the observed crumpling 
"transition" for fluid random surfaces is a lattice artifact in an asymptotically free theory. In addition 
the results here are another example of a simulation in the region beyond the C = 1 barrier of the 
continuum formulation which have found no obvious pathologies in either the graphs or the spin model. 
From the algorithmic point of view a comparison of the metropolis algorithm used here and a Wolff 
algorithm simulation revealed that there were unexpected equilibration problems on dynamical lattices 
for the Wolff algorithm as f3 was increased. In view of the results of |s) it would be most interesting 
to carry out a careful measurement of scaling on a dynamical lattice close to f3 = (where the Wolff 
algorithm is still well behaved) in order to see how the renormalization group flow compared with the 
fixed lattice model, and we are currently doing this. 

This work was supported in part by NATO collaborative research grant CRG9I009I. CFB is supported 
by DOE under contract DE-FG02-9IER40672 and by NSF Grand Challenge Applications Group Grant 
ASC-92I7394 DAJ is supported at LPTHE by an EEC HCM fellowship, EEC HCM network grant and 
an Alliance grant. The simulations were carried out on workstations at Heriot-Watt University and the 
Cray YMP-EL at LPTHE, Orsay. 
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Figure Captions 

Fig. 1. The specific heat for the various graph sizes simulated using a metropolis algorithm for the spins. 

Fig. 2. The flip acceptances, AL AF, and the probabilities of rings of length three, PR3, for a graph 
with N = 5000. 



■5 



